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In this Letter we illustrate the possible cyclic fermion pairing states across Feshbach resonances 
in optical lattices. In cyclic fermion pairing, the pairing amplitude exhibits an oscillatory behavior 
as the detuning varies. We estimate the quasi-particle gaps in different regimes of the resonances. 



Feshbach resonances in optical lattices were recently 
investigated in experiments [lj, |2( , after the earlier obser- 
vation of Feshbach resonances in traps 0, 0, 0- El- The 
fascinating phenomenon of lattice Feshbach resonances 
(LFRs) has attracted a lot of attention; both many-body 
physics and few-body aspects across LFRs have been 
studiedQ, H H G3 

[ill . Across LFRs fermionic atoms 
were observed to form shallow molecules, and superflu- 
ids should be present in optical lattices when tunnelling 
is permitted [2j. However, superfluids in LFRs can ex- 
hibit distinct properties. For instance, repulsive inter- 
actions between molecules 0, which have usually been 
overlooked in standard theoretical studies of superfluids 
near Feshbach resonances in traps, play a significant role 
in LFRs; the suprefluid long-range-order can be substan- 
tially renormalized by these interactions 

In this Letter, we study the pairing between fermions 
near LFRs, taking into account repulsive interactions be- 
tween bosonic molecules in the closed channel. We show 
that in the limit of large repulsive interactions between 
molecules, or small molecule band width, LFRs should 
exhibit a cyclic behavior as the detuning varies. Namely, 
the fermion quasi-particle energy gaps and critical tem- 
peratures of superfluids across LFRs exhibit an oscilla- 
tory behavior; the gap reaches maxima in a hard-core 
boson limit. 

Consider the following Fermi-Bose lattice model 
H = H f + H b + H bf ; 

H f = - t f,Wf(fLi<rfl*' + h - c -) 

<kl> ,rf,rf' ,<t 

H b = -t b ( b l b i + h - c -) 

<kl> 

+ y^(2ti - 2jj)b\b k + U b hbkjnbk - 1); 

k k 

H b f = ~Jo^2gw( b lfkrilfk v 'i+h.c.) (1) 

k , T) 

Here k, rj and a label lattice sites, on-site orbitals and 
spins; 77 = 1,2, ...M, a =j,|. fl w (f krj<T ) is the cre- 
ation(annihilation) operator of an atom (or fermion) at 



site k, with on-site orbital energy e v and spin a. b\{b]^) 
is the creation (annihilation) operator of a molecule (or 
boson) at site k (For simplicity, We assume there is only 
one bosonic orbital degree of freedom at each site). The 
boson number operator is n b k — b\ bk - tf,vv' an< ^ ^ are 
hopping integrals of fermions and bosons respectively and 
hopping occurs over neighboring sites labeled as < kl >. 
The single particle tunneling of bosons in optical lattices 
is assumed to be small, so that tb is much less than Ub, the 
on-site repulsive interaction between bosonic molecules; 
tf is much larger than the energy difference between two 
orbitals. Hb is the Hamiltonian previously introduced for 
atoms in optical lattices 0]; ^ has been used to study 
Mott insulator-supcrfiuid phase transitions which were 
recently observed in optical lattices [lij. 

Finally in narrow resonances, or in the large-N limit 
treated before [ij H^|. the matrix g vv > is proportional to 
S^i. However, in other limits the (/-matrix generally has 
off-diagonal terms 0] . Here we assume g vv ' to be an arbi- 
trary symmetric matrix so that the spatial wavefunction 
of pairs is symmetric and spin wavefunction antisymmet- 
ric. 

The self-consistent chemical potential is determined by 
the following formula, 

N = 2N b (2n - 2v; t b , U b ; A b (7o,t/,/x)) 
+N f (n;t f ;X f (j ,tb,Ub)) (2) 

where N, Nb and Nf are, respectively, the total number 
of fermions, the number of molecules or bound states of 
atoms and the number of atoms per lattice site. ^f( b ) 
is the effective interaction constant between fermionic 
atoms (molecules) mediated by molecules (fermions). We 
are interested in the limit where the coupling strength 70 
is small compared with U b so that in the calculation of 
N b (Nf) we can neglect the Ab(A/)-dependence. 

When t b is much smaller than U b , molecules can be 
considered to be in a Mott insulating ground state to 
the zeroth order in 70 Therefore N b is a step- like 

function of the detuning v. In the limit of zero t b , Nb is 
equal to an integer n whenever \2fi — 2v/Ub — L\ < 1/2; 
Nb has two degenerate solutions when 2fj, — 2v/Ub = I + 
1/2 and the value of N b jumps from / to / + 1 at that 
particular detuning signifying a discontinuous transition 
between two plateaus in the N b versus v curve. 
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FIG. 1: (Color online) Schematics of the occupation of atoms 
and molecules near a Feshbach resonance when Ub is large, 
(a) is when X is not inside the window of [0, 1]; (b) is for 
X £ [0,1]. (X is defined in Eq.(3).) In both (a) and (b), 
fully shaded lobes stand for completely filled molecular Mott 
bands and empty lobes for unfilled Mott bands; the partially 
filled lobe is for a partially filled band. The patterned area to 
the right of the vertical axis in (a) or (b) is the partially filled 
atomic band. The distance from the bottom of the atomic 
band and the bottom of the lowest molecule Mott band is the 
detuning \v\. 





(b) 



A finite but small value of t b won't change the gen- 
eral structure of the step-like function, except that the 
transition from the ith plateau to the (/ + l)th plateau 
becomes smooth. Generally, we have 



N b (2fj,-2v;t b ,U b ) 



Int 



2fi-2v f2fj,-2v-Int(2fi-v/U b )U b 1 



U b 



F 



1. 



F(X) = { f(X), 
0, 



2zt b 

< X < 1; 
< X < 0. 

4zt b — — 



(3) 



Here IntY = I if -1/2 < Y - 1 < 1/2; X is the distance 
(in energy space as shown in Fig^l from the chemical 
potential to the bottom of the nearest Mott band mea- 
sured in the unit of 2zt b (where z is the coordination 
number). The filling factor f(X) is a continuous func- 
tion of X varying from at X = to one at X = 1. 
Its detailed form depends on the band structure, which 
is not specified here. 

Taking into account this general property of N b , one 
easily finds that the solution to Eq.© is also a step-like 
function of v (see FigEJ ■ Whenever the detuning is such 



FIG. 2: (Color online) Schematics of solutions obtained in 
this letter for N/2 = 3 + e (0 < e < 1). (a) the ^-dependence 
of the chemical potential \i for atoms. The thick line is for 
the limit when Ub is much larger than tb\ the dashed-dotted 
line represents the limit when tb approaches zero. The thin 
blue line is the typical chemical potential in the weakly inter- 
acting limit when Ub is much smaller than tb- (b) The thick 
curves show the v-dependence of the BCS quasi-particle gap 
when Ub is large. The gap saturates at the value given after 
Eq. llH when approaching the dashed lines. Dashlines locate 
the position of the center of a Mott band along the ii-axis. In 
both (a) and (b), the thin blue curves are for the limit when 
Ub vanishes. 



that - j^<X<0orl<X<^, varying the detun- 
ing v does not change N b or Nf and therefore the chem- 
ical potential for fermions must be a constant 17]. This 
leads to plateaus in the self-consistent chemical poten- 
tial. The first plateau represents the limit of extremely 
weakly interacting superfluids, i.e. the detuning is posi- 
tive and much larger than the fermion energy of atoms. 
The chemical potential in this limit barely depends on 
the detuning v or coupling constant 70; its value is ap- 
proximately equal to the Fermi energy of a noninteract- 
ing Fermi gas, i.e. /ii ~ ep ; ep is the Fermi energy of 
a gas with N atoms per site. The value of the chemical 
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potential on the ith plateau is given as 

M? /a = /i? /8 -(i-l)(4) 3/ V< W (4) 

where e Q F is the Fermi energy of an atomic gas with two 
atoms per lattice site. The number of plateaus in the 
chemical potential i max is determined by the density of 
atoms; for N/2 — n + e (n is an integer and < e < 1), 
the value of imax and the number of valleys in Fig|2fb) 
are n + 1 . 

The transition between two plateaus takes place when 
the self-consistent chemical potential is located inside a 
Mott band; i.e. when X is between zero and one. The 
width of the transition region when the detuning varies 
is proportional to the band width 2zt}>- The distance 
between two consecutive transition regions is Vi+i — Vi = 
Ub', Vi is the detuning at the ith transition region when 
X(Vi) = 1/2. 

The trial wavefunction for the fermion sector is 

\g. s . >= nKq)+^(q)/^)/a(-* >; ( 5 ) 

q 

q is the crystal quasi-momentum and u 2 + X)r;{ v r)£ = ^" 
The solutions when chemical potentials are located on a 
plateau are distinctly different from the solutions when 
the chemical potentials are in a smooth transition re- 
gion (Fig|2a)). We will discuss these two limits sepa- 
rately. When X (defined in Eq.(3)) is either larger than 
unity or smaller than zero, the chemical potential falls 
on a plateau. Unoccupied molecule states, with which 
fermions in extended states have resonances, are well 
above the self-consistent chemical potential. The dif- 
ference between the energy of the lowest molecule state 
and the chemical potential depends on the detuning v 
or X. The energy difference 28\{X)ztb{> 0) is equal to 
-2Xztb if > X > -Ub/iztb and equal to Ub - 2Xzt b if 
1 < X < U~b/4:zt b ; typically it is a fraction of Ub- When 
S2(X) ^> 1, one obtains the effective on-site attractive 
interaction between atoms mediated by molecules 



v *'** = - x t**M> x * = 2MxjzT b ' (6) 

Minimization of the free energy < Ti > — fiN leads to the 
following equations for (u,v ri ^)(V is the volume): 

2uv^ v > v (q) = A(g^u 2 - g^^v^^v^); 

A = y ftrf«*(q)«(q)- ( ? ) 

Here &,„'(q) = e w (q) 

— Mi e »)i)' i s the kinetic energy 
tensor. The lattice constant has been set to unity. 



Solutions to this equation depend on the details of g V £ 
and the dispersion ^ vr) i (q) . In the limit when the different 
bands are degenerate, i.e. £ w '(q) = <W£(q), one obtains 
= gri^v; (u, v) are the usual BCS solutions, 

2uv _ A 
u 2 — v 2 £(q) ' 

Here E(q) — -\/C 2 (q) + |A| 2 is the quasi-particle energy; 
E w '.9r^' = L The energy gap A - i / exp(-2/A / ) is 
exponentially small. A/ is a dimcnsionless interaction 

constant estimated to be C 2 s 1 (x)zt b t f ' we nave assumed 
the atomic density of states at the chemical potential is 

1 /3 

cNy/tf, c is a quantity of order unity which depends 
on the details of the band structure. 

When < X < 1, the chemical potential is between 
two plateaus and varies continuously as a function of 
X within this window. In the mean field approxima- 
tion one obtains the following modified self-consistent 
equation, A = 70 < b >, -5 2 (X)2zt b + A < b > 2 = 
7o X) q (l/2^(q) — l/2e(q)); X(U b ) is an interaction con- 
stant appearing in an effective theory and is of order of Ub 
[13; 6 2 {X) is equal to X when < X < 1. In the lead- 
ing order, one neglects the 70-dependence of < b > and 
A is linearly proportional to 70. So in the dilute limit 
(X < 1), < b > is of order of y / S 2 (X)2zt b /U b . How- 
ever once the chemical potential of interacting molecules 
X(Ub)Nb becomes comparable to the Fermi energy of a 
Fermi gas with the corresponding density (which scales 
as N^ S tb), the mean field approach fails; this sets the 
limit on the number of atoms per lattice site, which has 
to be much less than (t b /U b ) 3 . The window where mean 
field theory is applicable thus shrinks to zero as t b /U b 
approaches zero. 

To get an estimate of < b > in the large- U b limit, it 
is useful to map the model to an XXZ model. Since the 
completely filled Mott bands are irrelevant to the discus- 
sions of superfluids, we only need to take into account the 
partially filled band. When tb is much less than Ub, the 
molecules can be treated as hard-core bosons at energy 
scales much lower than U b . For hard-core bosons, the on- 
site Hilbert space is one boson and no bosons. The aver- 
age number of bosons at each lattice site varies from zero 
to one depending on the position of chemical potentials. 
This suggests an effective pseudo spin S — 1/2 subspace 
(on-site): |cr^ = -1 >= b\\vac > b , \a% = 1 >= \vac > b . 
\vac >b is the vacuum of bosons. Furthermore, in this 
truncated space one can easily verify that three Pauli 
matrices can be represented as o~t— > 1 — 2b\pk, o~~£ — > bk 
and <j- -> b\; a ± = {a x ± ia y )/2^ 0. 

A state with cr| = 1 corresponds to a vacuum of 
molecules, and er| = — 1 a Mott state with one molecule 
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per lattice site. On the other hand, in translationally in- 
variant states where <y% has a finite expectation value in 
the XY plane, the condensation amplitude of molecules 
is nonzero; i.e. 

< b k >=< °l > ■ (9) 

So any state with all < > pointing along a specified 
direction that differs from the z-axis breaks the in-plane 
(XY type) rotational symmetry of the pseduo spin and 
corresponds to a superfluid. The effective Hamiltonian 
for molecules can then be written as 

H eff = "I E + °W + K E *i (10) 

<kl> k 

Here h\ = (X — l/2)2zi&. The Hamiltonian is invariant 
under a rotation around the z-axis or has an XY sym- 
metry which corresponds to the usual U(l) symmetry of 
a boson field. 

One easily finds the ground state of the above Hamil- 
tonian using a semiclassical approximation^^ and gains 
insight into superfluids in the large-L^ limit. In the 
semiclassical approximation, one assumes that the ori- 
entation of a pseudo spin is specified by a polar angle 

and an azimuthal angle <& as: < a x >— sinOcos$, 
< a y >= sin 9 sin $ and < a z >— cos 6. 

Minimizing the mean field energy with respect to O 
leads to the following solution for O. O remains at as 
X becomes negative, corresponding to an empty Mott 
band. As X becomes larger than 1, O is equal to it 
representing a Mott state with one molecule per site. For 

1 > X > 0, one obtains cosO = 1 — 2X, or sinO = 
2y/X(l — X); note that in this limit all solutions with 
a given O but arbitrary $ are degenerate and break the 
usual U(l) symmetry of a superfluid. Following Eq.(9) 
and Eq.(10), one arrives at A = 70 sin 9 exp(— i$), or in 
terms of X 



A = 7o 2 y/X(l - X) exp(-£$). (11) 

At X = 1/2 or half filling, sin 9 = 1 and the modulus 
of the expectation value < b > (proportional to sin 9) 
reaches its maximum. This leads to the maximal value 
of the BCS gap, |A| = 70. 

The calculation of T c in the large J7& limit is tricky 
because of large quantum fluctuations present when the 
chemical potential falls in the middle of a Mott band. 
Here we estimate T c again in two limits: a) when the 
chemical potential falls on a plateau corresponding to 
various valleys in Fig|21 (b); b) when the chemical po- 
tential is between two plateaus. In case a), fermions 
are weakly interacting and we expect the BCS theory 
is valid. T c scales as A, the solution to Eq.(8) and is 



exponentially small when 70 is much smaller compared 
to Ub- In case b), T c should be approximately the BEC 
transition temperature of molecules in a partially filled 
Mott band, which is proportional to e 2 ^ 3 i/ (e(< 1) is the 
filling fraction.). This temperature differs from the quasi- 
particle gap calculated in Eq. ifTTt and is also independent 
of the coupling constant 70. It reaches a maximal value 
when X = 1/2 because of the particle-hole symmetry. 
Therefore, we expect an oscillatory behavior in the crit- 
ical temperature as well. At a temperature lower than 
the BEC temperature but higher than the exponentially 
small BCS transition temperature when the chemical po- 
tential is between two Mott bands, the ground state is 
expected to alternate between a fermionic superfluid and 
a normal fluid as the detuning is varied. This might be 
observed in future experiments on LFRs. I would like to 
thank T. Leggett, K. Madison and G. Warner for discus- 
sions. This work is supported by a grant from UBC, a 
Discovery grant from NSERC, Canada and by the A. P. 
Sloan foundation. 
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